Dark Matter Detection with Polarized Detectors 
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We consider the prospects to use polarized dark-matter detectors to discriminate between various 
dark-matter models. If WIMPs are Dirac particles, with a particle-antiparticle asymmetry, and 
if those particles have neutrino-like parity-violating interactions with ordinary matter, then the 
recoil-direction and recoil-energy distributions of nuclei in detectors will depend on the orientation 
of the initial nuclear spin with respect to the velocity of the detector through the Galactic halo. 
If, however, WIMPS are scalars, Majorana fermions, or are matter-antimatter symmetric, no such 
signal can arise. Since the amplitude of this polarization modulation is fixed by the detector speed 
through the halo, in units of the speed of light, exposures several times larger than those of current 
experiments will be required to be probe this effect. 

PACS numbers: 98.80.-l5 



Although dark matter has been known for several 
decades to dominate the mass budget of galaxies, its 
particle nature is still mysterious. The coincidence be- 
tween the interaction strength required for an early- 
Universe rehc to have the right cosmological density and 
the electroweak interaction strength motivates the idea 
that dark matter is composed of some weakly-interacting 
massive particle (WIMP) [1-4]. However, WIMPs con- 
stitute a broad class of dark-matter candidates, includ- 
ing heavy fourth-generation neutrinos, various supersym- 
metric particles, particles in models with universal ex- 
tra dimensions [5], etc.; they may be scalar particles or 
fermions, and if fermions, Majorana or Dirac particles. 
The precise nature of the couplings of dark matter to 
ordinary particles varies considerably among the models. 

An array of searches for WIMPs is now underway, but 
terrestrial direct-detection experiments, designed to de- 
tect nuclear recoils from collisions with dark-matter par- 
ticles in the Galactic halo, provide likely our best hope 
to detect dark matter [6]. These detectors measure the 
energy of the nuclear recoils; the spectrum of such re- 
coil energies can then be used to discriminate a WIMP 
signal from background, and in case of detection, to con- 
strain WIMP parameters and discriminate between dif- 
ferent WIMP candidates. It has also been suggested [7- 
9] that the direction of the nuclear recoil can addition- 
ally be used to distinguish backgrounds and to constrain 
dark-matter parameters, and this approach is now being 
implemented experimentally [10]. 

However, there is yet another handle these experiments 
can exploit: the spin polarization of the detector nuclei. 
If WIMPs arc scalar particles, then their interaction rate 
cannot depend on the polarization of the nuclear spins. 
Likewise, Majorana dark-matter particles will not be po- 
larized, and so their interactions with a detector cannot 
depend on the detector polarization. Suppose, however, 
that dark-matter particles are Dirac fermions and that 
there is a particle-antiparticle asymmetry [11]. If so, and 



if these particles have a parity- violating interaction with 
ordinary matter, then the total detection rate as well as 
the recoil and energy /direction distribution can depend 
on the polarization of the target nuclei. Thus, measuring 
the polarization dependence of these distributions may 
help discriminate between backgrounds and shed light 
on the WIMP's particle nature in case of detection. 

To illustrate, suppose a dark-matter particle x of mass 
TO^ interacts with a nucleus N of mass mjv via a four- 
Fermi interaction of the form, 



(1) 



where G is a dimension-2 coupling constant, and a, b, c, 
and d are real parameters. This interaction gives rise to 
a differential cross section, 

^ = A + Biv-s) + B'{v'-s) + 0{v^), (2) 
dh 

for a WIMP of incident velocity v to scatter a nucleus, 
initially at rest with polarization s, to a recoil energy E 
and a final WIMP velocity v'. Here, 
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FIG. 1: An incident dark-matter particle with initial velocity 
V scatters from a detector nucleus with polarization s and 
recoil momentum q. Earth's velocity Ve is chosen to lie along 
the 2:-axis and the recoil rate in Eq. (7) is differential in both 
polar angle a and azimuthal angle /3. 
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The cross section for antiparticles (as opposed to parti- 
cles) is obtained (assuming CP is conserved) by replacing 
B — )■ —B and B' —J- ~B'. The dependence of the cross 
section, Eq. (2), on the dot product of a polar vector 
(v or v') with an axial vector (s ) is a manifestation of 
parity violation. Parity violation requires that at least 
three of the parameters a, b, c, and d be nonvanishing. 
While we leave a discussion of detailed models to fu- 
ture work, we do note that currently acceptable versions 
of Dirac-neutrino dark matter [12, 13] have such cou- 
plings, and such couplings may also be found in recent 
models of parity-violating composite dark matter [14]. 
(Recall, however, that the polarization-dependent effect 
we discuss below requires a particle-antiparticlc asym- 
metry, an ingredient not necessarily employed in these 
particular papers.) In the maximal-parity-violating case, 
a ^ -b = c = -d = 1/2, A = B = G^mNiSnv^)-^ and 
B' = 0. We will assume this case for our numerical work 
below. 

We now calculate the distribution of recoil energies 
and directions assuming that the detector moves with 
velocity Ve (which throughout we will take to be along 
the z axis) through the Galactic halo. If we were to ig- 
nore directional information, then we would simply cal- 
culate a (single) differential event rate dR/dE. If we 
were to consider the combined recoil energy/direction 
distribution for an unpolarized detector (as considered in 
Refs. [7, 8]), then we would calculate a double-differential 
rate dR/dE/d cos a, where cosa = Ve- and q is the di- 
rection of the nuclear recoil. If, however, the detector 
has a spin polarization s, which we take to be in the x-z 
plane, at an angle d from the z axis, then there may be an 
additional dependence on the azimuthal angle /3, about 
the z axis, between the recoil direction q and s. We must 
therefore in this case calculate a triple-differential event 
rate dR/dE/dcosa/d(3. 



FIG. 2: Contour plot of the differential event rate, Eq. (12), in 
units of events/kg/day/keV/sr as a function of cosa, where 
a is the polar angle, and the azimuthal angle /3, for fixed 
recoil energy E = 30 keV. We take the angle between the 
polarization and the detector velocity through the Galactic 
halo to be i9 = 90°. We also take G = (100GeV)"^ and 



= 100 GeV. 



We begin by writing the triple differential cross section 
dR/dEdfl as by demanding v ■ q = q/2fiv, where fi is the 
reduced mass of the WIMP-nucleus system, and q is the 
recoil momentum. Then, 
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S ( COS7 — 7^-^ ) = — 37;<5(w • q — Vq), 

(6) 

where Vq ~ q/'2fJ- is the nuclear recoil velocity, dfl = 
d cos ad/3 is a differential solid angle, and 6(x) is the 
Dirac delta function. 

The differential event rate [8] per unit detector mass is 
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where is the local WIMP number density, and f{v) is 
the dark-matter velocity distribution in the lab frame. To 
isolate the polarization dependence in the general case, 
we subtract signals with opposite spin orientations. 
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We assume a standard Maxwellian halo, for which 
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where vq = 200 km/s is the halo velocity dispersion and 
Vesc = 500 km/s is the escape speed from the Galactic 
halo. 

To evaluate the integral in Eq. (8), we choose the 
spin vector to be in the direction s = (sini9, 0, cosi?), 
the initial WIMP-velocity direction to be -C = 
(sin6' cos (/), sin6'sin0, cos 0), and the recoil direction to 
be g = (sin a cos /3, sin a sin /3, cos a) . We then use the re- 
lation S{g{v)) = S{v — vi)/\g'{vi)\ for the Dirac delta 
function, where g'{v) denotes differentiation with re- 
spect to V, and vi satisfies g{vi) = 0. This occurs 
when V = Vq/g{9,(f>,a, 13) = ui, where g{9,(j)^a, (3) ~ 
cos a cos 9 + sin a sin 9 cos{(j) — /3) . There is cos 9 depen- 
dence in the integrand in Eq. (8) through the cos 9 depen- 
dence of f{v), and there is 9 and (/) dependence through 

v' ■ s = v's [cos cos I? + sin sin I? cos ^] . (11) 

Performing the v integral in Eq. (8), we obtain 

dAR _ G^n^ f^^ v'^ v' ■ s 

dEdcosadl3 ~ 8A^7r2 ^ 15(6*, </>, a, /3)| 

xe(u'-lWn)e(«esc-t''), (12) 

where v' and g both carry dependence on 9 and 0. Here, 
Vmin = {rnNE/2ii^y/^ is the minimum WIMP velocity^ 
required to produce a nuclear recoil energy E, and <d{x) is 
the unit step function. The polarization-induced break- 
ing of the azimuthal symmetry about Ve prevents us from 
going further analytically, as can be done otherwise [8], 
but the remaining double integral is straightforward to 
evaluate numerically. 

Figs. 2, 3, and 4 show numerical results for 
dAR/ dE /d cos a/ d(3, the polarization-dependent part of 
the triple-differential event rate, as a function of the 
recoil energy i?, polar angle a, and azimuthal an- 
gle f3 of the recoil nucleus. The key feature here is 
the (3 (azimuthal-angle) dependence in Figs. 2 and 4 
which would not arise without parity-violating inter- 
actions and particle-antiparticle asymmetry. The de- 
pendence of dR/dE/d cos a also depends on the polar- 
ization, as shown in Fig. 3, even without azimuthal- 
angle information. Fig. 5 shows the differential re- 
coil spectra [1/ R){dR/ dE) (blue curve, color online) and 
{1 /\AR\){d\AR\/ dE) (red curve) obtained by integrating 
over the angular dependence in Eqs. (7) and (8) respec- 
tively. Thus, a polarization-dependent detection rate can 
be sought even without directional information. 

Since the detector will be fixed in the Earth frame, the 
daily revolution of the Earth provides a natural mod- 
ulation of the detector orientation with respect to the 



The analysis can be generalized to inelastic dark matter [15] by 
using Ui^in = [A -I- (niME/n)] {2mt^E)-^/^. 
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FIG. 3: Contour plot of the differential event rate, Eq. (12), 
in units of events/kg/day/keV/sr as a function of cos a and 
the recoil energy E, for fixed azimuthal angle /3 = 0. We 
assume = 90°, G = (100GeV)-^ and = 100 GeV. 
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FIG. 4: Contour plot of the differential event rate, Eq. (12), 
in units of events/kg/day/keV/sr as a function of E and /3 
for fixed polar angle a = 45°. We assume •& = 90°, G = 
(100 GeV)"^ and = 100 GeV. 



Earth's velocity through the halo. This can be used 
to isolate the dark-matter signal from systematic exper- 
imental effects (e.g., variations in sensitivity of the dete- 
tor with recoil direction) that might mimic such a signal. 
The rotation of the Earth around the Sun provides an 
additional systematic check. 
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FIG. 5: Plot of the differential event rate {1/ R){dR/ dE) 
(red, color online) and polarization-dependent differential 
event rate f /(| Ai?|)(i| Aii| /dE (blue) after numerically in- 
tegrating over a and jS. As before, we assume i? = 90°, 
G = (100GeV)~^ and = 100 GeV. 

A quick estimate shows that the experimental exposure 
needed to isolate polarization dependence may be feasi- 
ble. Consider, for example, a WIMP candidate with to- 
tal detection rate i?tot • Since the polarization modulated 
amplitude is velocity suppressed, AR/Rtot ~ w ~ 10~^, 
it is necessary to observe approximately 3 x 10^ total 
events for a 3ct discovery. Thus, if the DAM A Nal re- 
sults are due to dark matter with a total scattering rate 
^tot 0-5 kg • day with a cumulative exposure of ^ 1 
ton-yr ~ 4 x 10^ kg • day [16], then a future experi- 
ment will need to amass an exposure roughly an order 
of magnitude larger than DAMA's to observe this sig- 
nal. While larger than those of current detectors, this 
exposure is within the scale of those considered for de- 
velopment within the next decade [17]. 

While we have focussed here on WIMP-nucleon scat- 
tering, similar ideas may apply to detection schemes 
based on dark-matter-electron scattering [18], in which 
the signal may conceivably be enhanced by the larger 
scattering rates associated with more abundant lower- 
mass dark-matter candidates. The increased relative ve- 
locity between the dark-matter particle and the target 
particle (an atomic electron) may also enhance the po- 
larization dependence. 
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